Directed Chaotic Transport in Hamiltonian Ratchets 
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We present a comprehensive account of directed transport in one-dimensional Hamiltonian sys- 
tems with spatial and temporal periodicity. They can be considered as Hamiltonian ratchets in the 
sense that ensembles of particles can show directed ballistic transport in the absence of an average 
force. We discuss general conditions for such directed transport, like a mixed classical phase space, 
and elucidate a sum rule that relates the contributions of different phase-space components to trans- 
port with each other. We show that regular ratchet transport can be directed against an external 
potential gradient while chaotic ballistic transport is restricted to unbiased systems. For quantized 
Hamiltonian ratchets we study transport in terms of the evolution of wave packets and derive a 
semiclassical expression for the distribution of level velocities which encode the quantum transport 
in the Floquet band spectra. We discuss the role of dynamical tunneling between transporting is- 
lands and the chaotic sea and the breakdown of transport in quantum ratchets with broken spatial 
periodicity. 

PACS numbers: 05.60.-k, 05.45. Mt 



I. INTRODUCTION 

Hamiltonian systems with a mixed phase space remain 
a challenge within the field of nonlinear dynamics, both 
classical and quantum. This is usually attributed to the 
intricate, typically self-similar structure of phase space 
in these systems. There exist, however, more tangible 
effects which also require a coexistence of regular and 
chaotic dynamics but no particular fine structure. A 
prominent example is directed transport: An elementary 
yet decisive consequence of a mixed phase space is the 
existence of distinct regions which support qualitatively 
different dynamics and do not communicate with each 
other. Directed transport may arise locally in regular 
components of phase space. As a consequence of a global 
sum rule, and in the absence of certain symmetries, it 
can then be conferred to the chaotic component, as we 
will show in this paper. 

Chaotic transport in extended Hamiltonian systems is 
usually associated with undirected diffusion: The width 
of the spatial distribution Ax grows with time as some 
power law (Aa;)^ ^ with a between and 2. Only 
recently it has been discovered that even in the absence 
of a mean external gradient, chaotic diffusion in driven 
Hamiltonian systems can be accompanied by a directed 
drift. The corresponding ballistic component of trans- 
port [1] may surprise on first sight, since a hallmark of 
chaos is the decay of all correlations including an effective 
randomization of the velocity with time. However, this 
implies only that the mean velocity of a typical chaotic 
trajectory approaches an asymptotic value which is char- 
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acteristic of the chaotic phase-space region as a whole. In 
the absence of additional symmetries there is no general 
reason requiring this asymptotic mean velocity to be zero. 

In fact, as we shall argue in Sec. II, in systems with 
a mixed phase space a sum rule requires chaotic trans- 
port to compensate for the directed transport possibly 
occurring in regular phase-space regions [2, 3]. An im- 
portant conclusion (Sec. II E) is that the ballistic chaotic 
transport has nothing to do with internal structures of a 
chaotic phase-space component such as cantori or other 
partial transport barriers. All these complicated sub- 
structures, leading, e.g., to Levy walks and anomalous 
diffusion in Hamiltonian ratchets [1, 4, 5], need not be 
considered in detail in order to understand that ballistic 
transport dominates for long times. 

Deterministic ballistic transport due to a dynamical 
restriction of trajectories to certain phase-space regions 
has been observed before in dissipative systems [6] , where 
phase-space volume is contracting with time. This mech- 
anism is close to the concept of stochastic ratchets (Brow- 
nian motors), i.e., systems that generate directed motion 
from non-equilibrium noise [7-9]. The analogy suggests 
to speak of deterministic ratchets. 

Throughout this paper, we disregard dissipation. Its 
absence, however, not only renders it more difficult to 
achieve directed transport, since the natural time arrow 
determined by dissipation is lost and has to be replaced 
by other mechanisms breaking time-reversal invariance. 
It even becomes a subtle task to define a ratchet in a 
Hamiltonian framework in the first place. Trajectories 
can maintain a memory of their initial velocity for an 
infinite time. Therefore a precise definition of a Hamilto- 
nian ratchet is not completely straightforward: The mere 
fact that in unbiased systems directed transport can ex- 
ist and survive for infinite time is trivial, just take a free 
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particle with some non-zero initial velocity vq 7^ 0. In 
this sense every extended Hamiltonian system would be 
a ratchet. 

Due to velocity dispersion an ensemble of free particles 
will also spread ballistically, i.e., as fast as its center of 
mass is transported. On the other hand, as pointed out 
above, there exist Hamiltonian systems where transport 
is ballistic, but the spreading is not. They are character- 
ized by a locking of the average velocity to a (non-zero) 
value which does not depend on the precise initial condi- 
tions as long as they are restricted to some finite phase- 
space region. For the purpose of the present paper we 
regard this property as the definition of a Hamiltonian 
ratchet. 

Even with this restriction it is possible to construct 
cases one would qualify as trivial realizations of directed 
transport: In the integrable system sketched in Fig. 1, 
for example, transport appears to be achieved by a mere 
change of frame. For the sake of simplicity of the defini- 
tion we do not attempt to formally exclude such cases. 
In what follows, however, we concentrate on extended 
systems with a mixed phase space where one has to un- 
derstand the interplay between regular and chaotic trans- 
port. 

Unless symmetries of the driving potential prevent it 
[1], Hamiltonian ratchets as defined above lead without 
average force to a non-zero mean velocity of an ensemble 
of particles which were initially at rest. The same applies 
also to the ratchets described in [10], although there is 
no velocity locking and ensembles of particles do spread 
ballistically. These systems are based on a mechanism 
that is different from the models discussed in [1-5, 11] 
and the present paper, and we will not consider them 
here. 

For Hamiltonian ratchets under the influence of an av- 
erage force we show in Sec. II F that uphill regular trans- 
port is possible. In contrast, even an infinitesimal aver- 
age force destroys the chaotic drift and leads to downhill 
acceleration. 

It comes as a rather unexpected finding that Hamil- 
tonian ratchets have applications on macroscopic, even 
geophysical scales where apparently friction prevails [12]. 
Indeed, in hydrodynamics, even in the presence of dis- 
sipation, restricting the description to position space re- 
sults in a Hamiltonian form of the evolution equations if 
only the fiuid is incompressible. Specifically, in geophysi- 
cal applications, a periodic potential reflects the periodic 
boundary conditions on Earth with respect to longitude, 
while an asymmetry in the transverse coordinate is im- 
plied by the dependence of the Coriolis force on latitude. 

Going in the opposite direction, Hamiltonian ratch- 
ets are to find applications on scales where quantum ef- 
fects become important. For example, in semiconductor 
nanostructures employed to investigate solid-state ratch- 
ets [13, 14] such effects were observed. A Hamiltonian 
ratchet with negligible dissipation can be realized on this 
basis if the structure size is further decreased, such that 
electronic motion occurs in the ballistic regime. But this 



will even enhance quantum corrections. 

In Ref. [2] it was concluded that quantum Hamiltonian 
ratchets can work if classical and quantum system are 
both spatially periodic such that the quantum system has 
a band spectrum. Detailing our findings, we will show in 
Section III B that quantum transport relies on the semi- 
classical correspondence between the dynamics of wave 
packets and that of classical distributions in phase space: 
As long as a wave packet, started in the chaotic region 
of phase space, say, remains predominantly restricted to 
this region, it will be transported with the classical mean 
chaotic velocity. Such quantum-classical correspondence 
can be attributed to the existence of different types of 
bands in the spectrum, with eigenfunctions concentrat- 
ing semiclassically on different invariant sets of classical 
phase space. Since this mechanism crucially depends on 
classical phase-space structures, it cannot be captured 
using a single- (or few-) band picture. Therefore our re- 
sults are not at variance with the absence of transport 
demonstrated within such an approximation [11]. 

However, also in the semiclassical regime non-classical 
processes like tunneling are possible which allow transi- 
tions between invariant sets of classical phase space. In 
Section HI C we will address the question why this is 
compatible with quantum transport unlimited in time. 
Only when the exact periodicity of the quantum system 
is destroyed, the eigenfunctions governing the long-time 
dynamics ignore classical phase-space structures [15] such 
that ratchet transport becomes a transient phenomenon. 
We shall deal with this case in Section HID. 

In our conclusions (Sec. IV) we discuss in particular 
various ways of breaking the translation invariance of 
Hamiltonian ratchets and how this affects transport. 
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FIG. 1: (a) A trivial example for a Hamiltonian ratchet is 
a periodic potential which is jnoving at a constant velocity 
vb > such that V{x,t) — V{x — vst)- The system is in- 
tegrable, since it is time independent in the comoving refer- 
ence frame. Despite the conveyor-belt construction it is also 
unbiased, since the average force in a periodic potential is 
always zero, (b) shows the dependence of the asymptotic 
mean velocity iJ on the initial velocity vq under the assump- 
tion that the potential is non-zero only in negligibly small 
intervals. Particles with initial velocity close to vb, namely 
for m(vo — vsY /2 < Vmax, are trapped inside one well of the 
potential and have an asymptotic velocity v = v-g, independent 
of the precise initial conditions. 
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II. CLASSICAL HAMILTONIAN RATCHETS 

A. The Hamiltonian of the extended system 

We consider Hamiltonian systems in one dimension 
which are periodic and unbiased in the sense specified 
below. The Hamiltonian is of the form 

H{p,x,t) = Tip) + V{x,t), (1) 

where X and p are the canonically conjugate position and 
momentum and T{p) and V{x,t) denote kinetic and po- 
tential energy, respectively. 

We require that the dynamics be invariant under inte- 
ger translations of space or time and use dimensionless 
variables in which both periods are unity, i.e., we assume 
the following property: For any trajectory x{t) with ini- 
tial conditions x{ta) = .to, p{to) = Pa and any other 
trajectory x{t) with x(to + n) = xq + m, p{to + n) ^ po 
we have x{t + n) = x{t) + m for all t. 

In the simplest case this is realized by T{p) = p^/2 and 
a spatially and temporally periodic potential 

V{x,t + l) = V{x+l,t) = V{x,t), (2) 

but this is not a necessary condition: If the potential 
contains an additional term f{t) x we have V'{x+l,t) = 
V'{x,t) only, where V = dV/dx. Nevertheless, discrete 
translation invariance may be satisfied for the dynamics, 
see Section II F for an example. 

We shall refer to the system as unbiased, if the force 
—V averaged over space and time vanishes 

f dx [ dtV'{x,t)=0. (3) 
Jo Jo 

In Section II G we will also consider systems where the 
kinetic energy is a periodic function of p such as T{p) = 
cos27rp for electrons in a Bloch band. As we shall see, 
such systems are always unbiased. 



B. The phcise space of a unit cell 

Instead of the extended system represented by (1), the 
discrete translation invariance allows to consider an aux- 
iliary system restricted to a single unit cell by imposing 
periodic boundary conditions at x = 1, t = 1. Since in 
this paper both representations appear in parallel, we use 
different symbols ^ = x mod 1 and t = t mod 1 for the 
cyclic variables of the unit cell. 

It is a standard technique for driven systems [16] to 
treat time like a spatial coordinate such that a one- 
dimensional time-dependent system is mapped to a for- 
mally time-independent problem in two dimensions. For 
the unit cell the Hamiltonian obtained in this way is 

n{^,p,T,£)=T{p) + Vii,T) + £, (4) 



where £ is canonically conjugate to r. This ensures f = 
dH/d£ = 1. Since H is a conserved quantity, —£{t) 
can be interpreted as the energy AiJ which the system 
has gained from the driving up to time t. Moreover it 
becomes clear that the dynamics is restricted to a three- 
dimensional "energy shell", Ti = const., which is spanned 
by the variables ^, p and r (f is a function of these three 
variables and the constant W). 

The dimensionality can be reduced further by consider- 
ing Poincare surfaces of section at some constant r which 
eliminates the trivial flow in r-direction. In the following, 
we shall discuss the main features of such stroboscopic 
surfaces of section, relevant for transport in Hamilto- 
nian ratchets. For the moment we restrict the discussion 
to smooth potentials in the sense of the Kolmogorov- 
Arnol'd-Moser (KAM) theorem [17] and take as an ex- 
ample the Hamiltonian 




1 

FIG. 2: Typical stroboscopic Poincare section r = for a 
Hamiltonian ratchet with non-contractible KAM tori, main 
chaotic sea and regular islands. The lettered rectangular re- 
gions support initial distributions of particles for which the 
corresponding velocity distributions are shown in Fig. 3. 
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with 

Vo{x) = -^[sin(27ra;) + 0.3sin(47ra; + 0.4)] (6) 
0.(0 

and 

Vi(x) = --^[4.6 sin(27ri) + 2.76 sin(47rt + 0.7)1 . (7) 
5.76 

This corresponds to the parameter set (3) of Fig. 1 in 
Rcf. [1] when the spatial and the temporal period are 
scaled to unity. 

The stroboscopic Poincare section for this model is 
shown in Fig. 2. Wc can distinguish the following three 
types of motion, each corresponding to a characteristic 
signature in phase space and transport: 

(i) At high kinetic energies the ratchet potential can be 
considered a small perturbation acting on a free particle. 
For this intcgrable limit the trajectories are confined to 
invariant surfaces in phase space which have the topol- 
ogy of a torus. These tori are labeled by the conserved 
value of the momentum p and parameterized by the cyclic 
variables ^ and r. In the {^,p) plane of the stroboscopic 
Poincare section the tori would consequently appear as 
horizontal lines. 

The KAM theorem predicts the fate of a torus under 
a small perturbation. It depends on whether its winding 
number w is rational or not. The winding number is 
the ratio between the angular velocities along the two 
independent cyclic coordinates spanning the torus. In the 
present case, one of these coordinates is the time r and 
the corresponding angular velocity is unity by definition. 
For the other coordinate ^, the angular velocity on the 
torus is equal to the transport velocity in the extended 
system, measured in spatial unit cells per time period, so 
that w = v. 

Almost all tori have irrational winding numbers and, 
according to the KAM theorem, most of them survive an 
infinitesimal perturbation. This is visible in Fig. 2 at high 
\p\ where wc observe lines in the stroboscopic Poincare 
section which extend across the unit cell. Although the 
lines are deformed by the potential they represent intact 
tori of regular motion with irrational winding number 
(transport velocity). Motion proceeds on these tori in 
the initial direction, without turning points. As these 
tori cannot be continuously contracted to a point we will 
call them non-contractible. 

(ii) Tori with rational winding number 

w = u^/Ht, (8) 

Ux, Hi integer, which comprise a set of measure zero, 
are destroyed under an infinitesimal perturbation. De- 
tails of this effect are described by the Poincare-Birkhoff- 
Theorem [17]. Together with a small neighborhood, a ra- 
tional torus decays to a chaotic layer embedding new tori 
of regular motion. These tori have a different topology, 
however: They are contractible and appear as a set of /i^ 
regular islands in the stroboscopic Poincare section. The 



respective centers of the islands are formed by a single el- 
liptic periodic orbit with period /i^, i.e., this orbit has /Xj 
distinct intersections with the stroboscopic Poincare sec- 
tion before it starts repeating, shifted by Vx unit cells in 
the extended space. In Fig. 2, we therefore observe chains 
of regular islands which are sequentially traversed by a 
trajectory. The average velocity of the central periodic 
orbit and of all trajectories inside the island corresponds 
to the rational winding number TJ = w of the destroyed 
torus. If V ^ we speak of a transporting island. 

(iii) The chaotic regions surrounding the island chains 
at high \p\ are too small to be visible in Fig. 2. With in- 
creasing perturbation, however, the chaotic regions grow 
and may coalesce. In the vicinity of p = the ef- 
fective perturbation is strongest. As a result a large 
chaotic sea develops. With increasing resolution we find 
more and more islands embedded in this sea and more 
and more chains of transporting islands interrupting the 
strips where the intact KAM tori reside. Such islands 
need not be remnants of rational tori in the undrivcn 
system — they can appear and disappear at some finite 
value of the driving potential as a result of bifurcations 
of periodic orbits. Still, their transport velocity must also 
be given by a rational winding number. 

Conversely, we find more and more small chaotic re- 
gions located within the regular islands. Since they are 
confined to the phase-space region demarked by the out- 
ermost intact torus encircling the island, they share the 
same average velocity v = w, where w is the winding 
number of the island. 

The phase-space regions enumerated above are most 
adequately discussed in terms of invariant sets: a subset 
of phase space invariant as a whole under the dynamics, 
irrespective of any reshuffling possibly occurring inside. 
For example, any regular torus in the three-dimensional 
phase space of the unit cell is invariant under the dynam- 
ics. Trajectories initialized on the torus do not leave it 
and vice versa. This invariant two-dimensional surface 
separates the remaining phase space into two invariant 
sets of non-zero measure. Moreover, any region in phase 
space confined by a number of tori is an invariant set of 
the dynamics. In particular, this applies to the chaotic 
sea, which is bounded from below and above by two non- 
contractible KAM tori and by the outermost tori of the 
embedded regular islands. 

For our purpose the limitation of chaotic trajectories to 
a compact region of phase space will be crucial. In the ex- 
ample discussed above this is a consequence of the KAM 
scenario. In systems where the KAM theorem is not valid 
our theory applies as long as there is another mechanism 
leading to a compact chaotic phase-space component. An 
example of this type will be discussed in Section II G. 



C. Velocity distribution 

Although the system (4) is restricted to a single unit 
cell, it contains the complete information about trans- 
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port in the extended system (1). The velocity v = 
dH/dp = T' (jp) along a trajectory is the same in both 
cases provided the initial conditions are equivalent, i.e., 
^0 = xq mod 1 at f = 0. Therefore the velocity is the ap- 
propriate quantity to connect transport in the extended 
system to the unit cell and we describe transport in terms 
of the velocity distribution for an ensemble of particles. 
An ensemble is specified by a normalized initial distri- 
bution pq{S^{),pq,tq) in the phase-space unit cell. The 
variable Tq is part of the initial conditions since it mat- 
ters at which phase of the driving force a trajectory was 
started. It can indeed be physically meaningful to con- 
sider ensembles for which tq is not sharp, for example to 
model a situation where particles continuously enter the 
system. 

For any ensemble po and time t we define the time- 
averaged velocity distribution as 



PpoM 



r'+oo 



\ I dt' j dpo 

r Jo J-oo 



d^o / dro 
Jo 

)) (9) 



with the normalization J dv P{v) = 1. If we consider an 
ensemble in the extended system, initially localized at 
a; = 0, then at a later time t its spatial distribution will 
be given in terms of the velocity distribution by pt{x) = 
PpQ^t{x/t). For long times the center of mass moves 
with the mean velocity 



-t-oo 



dv V Po. 



(10) 



where the existence of Ppo,oo{v) = limt_>oo Ppo,t{v) is as- 
sumed. 
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FIG. 3: Distribution of time-averaged velocities, Eq. (9), for 
four different initial distributions a-d (see Fig. 2) . The chaotic 
distribution (a) was sampled by 10,000 trajectories, while for 
each of the regular distributions (b d) only 100 trajectories 
were used. For each trajectory the velocity was averaged up to 
t = 100 and the resulting distributions are displayed with solid 
lines (a-d). For (a), (d) also the distributions at t = 10,000 
are shown (bold lines). 



The behaviour of the velocity distribution is qualita- 
tively different for initial distributions pQ which are re- 
stricted to different invariant sets of the phase space. 
This is demonstrated in Fig. 3. We used as initial distri- 
butions the characteristic functions Xa,b,c,d of the rectan- 
gles marked in Fig. 2, approximated by a large number of 
trajectories with initial conditions distributed randomly 
inside the corresponding region. 

In the simplest case, po has support inside a regular 
island (distributions b, c in Fig. 3). According to the 
last section, the average velocity of all trajectories inside 
an island is equal to the winding number w of the island. 
Consequently we have 



Ppo,oo{v) =S{v-w) 



(11) 



and observe sharp peaks in Fig. 3b, c whose width is 
within the bin size of the histogram already at t = 100. 
Fig. 3c is an example for a transporting island, w; ^ 0. 
Any distribution po initialized inside this island will be 
transported ballistically with velocity w = ~1. At the 
same time the width of the distribution does not grow 
ballistically. As stated in the introduction, we consider 
this behaviour as the defining property of a Hamiltonian 
ratchet. 

For an ensemble initialized in the chaotic sea (Fig. 3a) 
the situation is similar. Although here the velocity dis- 
tribution shows an appreciable width at finite times, the 
comparison of t = 100 and t = 10,000 suggests that 
this width goes to zero as t — > oo. We can explain 
this behaviour using the concept of ergodicity. Ergod- 
icity means that for any function defined on phase space 
and for almost all trajectories the time average along the 
trajectory coincides with an average over the accessible 
phase space. It is usually assumed that this property ap- 
plies to the chaotic components of systems with a mixed 
phase space, although proofs of such a statement can be 
given only in exceptional situations [18]. For our pur- 
pose wc can use the velocity v = T'{p) as the function on 
phase space and obtain for any non-singular initial dis- 
tribution inside the chaotic sea, such as the rectangular 
region of Fig. 2a, 



(90, oo 



{y) =5{v- Dch) 



with the mean chaotic velocity 



Vo^ = K^ [ drd^dpT'ip). 

Jch 



(12) 



(13) 



The phase-space integral extends here over the whole 
chaotic sea of the spatio-temporal unit cell, and Vch = 
J^^ dr d^ dp denotes its volume. 
In the following section we shall discuss a method to 

evaluate Eq. (13). For the moment it suffices to say that, 
in the absence of specific symmetries, there is no gen- 
eral reason to expect that the chaotic velocity predicted 
by this equation is zero. Therefore, also the chaotic sea 
provides an example for Hamiltonian ratchet transport. 
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For both, regular islands and chaotic components the 
asymptotic velocity distribution is a (^-function which 
does not depend on the precise location of the initial 
phase-space distribution within the invariant set. The ve- 
locity distribution obtained from a region with surviving 
non-contractible KAM tori shows a fundamentally differ- 
ent behaviour, analogous to the case of a free particle: it 
maintains a finite width for t —>■ oo and a complicated 
internal structure (distribution d in Fig. 3). Moreover, 
the detailed properties of the asymptotic velocity distri- 
bution depend on the precise shape and location of the 
initial ensemble. Hence, according to our definition, non- 
contractible tori do not show ratchet-like transport. 



D. Transport for invariant sets and sum rule 

There is an interesting reformulation of Eq. (13) which 
allows to calculate the chaotic mean velocity in terms of 
regular trajectories only [2]. For any subset M of the 
unit cell, we define its contribution to transport, Tm, as 
phase-space volume times average velocity, 

Tm = Vmvm 

didpdTXMii,P,T)T'{p), (14) 



where Xm{£,,P,t) is the characteristic function of M. 
Note that in this definition M is not necessarily an in- 
variant set. However, if M denotes either the chaotic sea 
or a regular island, the phase-space averaged velocity vm 
can be identified with the asymptotic mean velocity of al- 
most all trajectories inside the invariant set, as described 
in the previous section. 

Transport has to be distinguished from the familiar 
concept of current which refers to the probability flow 
that passes per unit time through a surface dividing 
phase space. Here we are interested in transport along 
the x-direction. Therefore we consider the current at a 
point ^0- The value of the current depends on the posi- 
tion ^0 and the time tq where it is measured. In terms of 
the density Pt{^,p), it is given as 



/(6,ro) = 



+00 



dpPro{i0,p)T'{p). 



(15) 



In order to relate this current to the transport of an in- 
variant set M, Eq. (14), we substitute the density of the 
invariant measure 



Xm{^,P,t) 



A 



M-. 



(16) 



where Am denotes the area of M in a stroboscopic 
Poincare section. Integration of the density over one pe- 
riod of the driving leads to the time-averaged current of 
M at ^0, 

_ 1 f+oo 

Im{Co) = ^ dT dpxMiCo,P,r)T'{p) (17) 



where we have used the conservation of phase-space area 
in Hamiltonian systems, Am,t = Am- An additional 
integration over yields the relation between current in 
a;-direction and transport 



T 



M 



AmI 



(18) 



Here we have used that the time-averaged current is inde- 
pendent of the position ^o, as implied by the continuity 
equation for the invariant measure. Note that for this 
reason we could in principle define transport also with- 
out the ^-integration. 

By choosing the density as in Eq. (16) and weighting 
the contribution of each invariant set M by its area Am, 
we achieve that the resulting quantity, transport, is ad- 
ditive. Namely, with the definition (14), we have for the 
union of two or more disjoint sets, i.e., for M = [j-Mi, 
for all i ^ j. 



with Mi n Mj 



T 



M 



(19) 



We will apply this sum rule for transport to the layer 
in phase; spacx; which contains the chaotic sea and the 
embedded regular islands. It is bounded from below and 
above by two KAM tori. For simplicity we assume that 
they can be represented by two functions Pu/i{^,t). We 
find from Eq. (14) 

^ayer = / / dr / dpT'{p) 

Jo Jo ■'p\{i,T) 

= /'dC /'dr[T(p,(e,T))-T(pi(e,T))] 
Jo Jo 

= {T)n-{T)u (20) 

i.e., the transport of the layer is simply given by the ki- 
netic energy T, averaged over the two bounding KAM 
tori. In short, since the underlying phase-space distri- 
bution Xm(^,P, t) is fiat, the transport is determined by 
the outline defining the subset M. This applies to any 
subset of phase space confined by two non contractible 
tori. 

On the other hand, according to Eq. (19) the transport 
of the stochastic layer is equal to the contributions from 
the invariant manifolds it comprises 



{T)^-{T), = VM + J2ViVi 



(21) 



Eq. (21) can be used to predict the chaotic transport 
velocity. In practice this works as follows: 

(i) In the stroboscopic Poincare section we determine 
the location of the limiting KAM tori p^^i and the lo- 
cation of the limiting tori of all major regular islands i 
together with their winding munbers Wi. 

(ii) In order to determine the phase-space volumes en- 
tering Eq. (21) it is in fact sufficient to know the areas in 
the stroboscopic Poincare section: The Liouville theorem 
applied to the time-dependent Hamiltonian Eq. (1) [17] 
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ensures that such an area is conserved by the dynamics. 
The three-dimensional volume within the phase space of 
the unit cell is simply the area at any given moment in 
time, multiplied by the temporal period V = A x 1. Areas 
in the Poincare section are determined by approximating 
the corresponding invariant manifold by a polygon with 
corners obtained from running a trajectory on the outer- 
most torus. Numerically, an approximation to this torus 
can be found by zooming into the Poincare section. 

(iii) The kinetic-energy averages (T)u,i over the bound- 
ing KAM tori are obtained by sampling a torus with a 
long trajectory, and determining the integrals Eq. (20) 
numerically. Note that this is not equivalent to a time 
average over such a trajectory as the invariant density on 
the torus is not constant. 

(iv) Putting all the information together we find 



Nevertheless substructures inside the chaotic compo- 



Vch 



{T)^-{T),-j:^AiWi 



A, 



layer 



(22) 



Compared to the above procedure, the straightforward 
method of determining the chaotic transport velocity by 
running a very long trajectory has the disadvantage that 
its accuracy is hard to control. The trajectory must be 
long enough to sample the chaotic phase-space compo- 
nent ergodically, and there is no way to tell from a single 
trajectory whether this has been achieved with sufficient 
accuracy. The reason is that the chaotic component typ- 
ically contains partial barriers (cantori), which may ap- 
pear closed in a simulation over finite time. The error 
made by ignoring the phase-space region behind the par- 
tial barrier can in principle be arbitrarily large. Also 
the converse error is possible: For long simulations the 
accumulating numerical inaccuracy may drive a chaotic 
trajectory beyond an intact KAM torus. By using a stro- 
boscopic Poincare section such errors are substantially 
reduced. In the picture obtained from many relatively 
short trajectories, sampling the entire phase space, one 
can judge if there are two nearby chaotic regions which 
may actually form a single invariant set. It is then suf- 
ficient to increase the resolution selectively in a small 
portion of phase space, which is possible with relatively 
small computational effort. 



E. Chaotic transport and Levy walks 

Equation (22) shows that the basic mechanism under- 
lying chaotic ratchet transport is the existence of KAM 
tori and regular islands which prevent a chaotic trajec- 
tory from sampling the whole classical phase space. Un- 
less there are special symmetries, the velocity average 
over the chaotic sea is generically non-zero and it is de- 
termined solely by the boundaries of this invariant set. 
Besides ergodicity, no reference to any details of the dy- 
namics within the chaotic set is needed to explain and 
quantitatively predict the observed asymptotic chaotic 
transport velocity. 
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FIG. 4: (a) Stroboscopic Poincare section at r = for the 
system of Eq. (24). (b) For an initial distribution po ~ Xch 
the distribution of time-averaged velocities is shown at vari- 
ous times t. As t — > 00 it evolves to a narrow peak around 
the asymptotic mean velocity (dashed line for t = 10^). (c) 
From all distributions shown in (b) the average velocity v 
is computed after the contributions from the balhstic chan- 
nels have been removed by restricting P{v) to the interval 
—28 < V < +18. The resulting values (dots) are for all times 
close to the asymptotic mean velocity (dashed). 
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nent of phase space in general do exist and leave their 
hallmark in transport properties. Levy walks, in par- 
ticular, have attracted some attention in the context of 
Hamiltonian ratchets [1, 4, 5]. These are the episodes 
when a chaotic trajectory is trapped in the vicinity of 
a transporting island, close to the hierarchical struc- 
ture of smaller and smaller islands and surrounding can- 
tori. Such hierarchical regions are virtually unavoidable 
in a mixed phase space (for remarkable exceptions see 
[18, 19]). In the context of ratchets they were termed 
"ballistic channels" [4, 5] and are frequently located in 
the vicinity of the KAM tori confining the chaotic sea 
from below and above, i.e., in regions of relatively high 
velocity. Therefore Levy walks are easily observed in nu- 
merical transport experiments. Some care must be taken 
to avoid the wrong conclusion that ballistic channels and 
Levy walks are necessary for the existence of substantial 
chaotic transport or can completely account for it. 

To study this question in some detail, let us start 
from the sum rule Eq. (19) and decompose the chaotic 
transport into contributions from disjunct subsets of the 
chaotic sea C = \Jj Cj. We have Vchi^ch = J2j Vi^j 
Vch = J2j Vi such that 



(23) 



Because of ergodicity inside the chaotic component the 

phase-space volumes Vj in Eq. (23) can be replaced by 
the fraction of time a typical chaotic trajectory spends 
inside subset j or, equivalently, by the probability to en- 
ter subset j times the average survival time in it. Doing 
so we immediately arrive at a formula similar in spirit to 
Eq. (3) of Ref. [4] or Eq. (6) of Ref. [5]. At the same time 
it is still exact and does not depend on the character of 
the subsets j used to subdivide the chaotic region. As in 
Rcfs. [4, 5], this decomposition can, e.g., consist of a few 
prominent ballistic channels and some remaining chaotic 
"bulk" region. Our main point is here that in general it 
is not possible to approximate this remainder by an undi- 
rected and purely diffusive dynamics, i.e., to set Vj = 
for the corresponding subset in Eq. (23). 

For this purpose we will follow the analysis suggested 
in Refs. [4, 5] but apply it to a model with different pa- 
rameter values. The Hamiltonian is 



H{p,x,t) 



P 



271 cos(27ra;) 



(24) 



+{2Trf X 2 cos(27rt) - 4 cos (^Trt -h ^) 



and the stroboscopic Poincarc section (Fig. 4a) shows 
the typical features discussed in Section II B. The ve- 
locity distribution of the chaotic component is shown in 
Fig. 4b for various times. In contrast to Fig. 3c we have 
chosen here an ensemble of initial conditions po ~ Xch 
uniformly covering the entire chaotic sea. Numerically 
this has been achieved relying on ergodicity: We run 
a single long chaotic trajectory x{t) (0 < i < 4 • 10^) 



and used x{t') with t' = 0, 1, 2, . . . as the initial condi- 
tions of the ensemble. For each such initial condition 
Vt = [x{t' + t) — x{t')]/t is the velocity averaged over 
a time span t. For fixed t the probability distribution 
P{vt) is shown in Fig. 4b. It is equivalent to the propa- 
gator used in Ref. [5] for visualizing internal details of the 
chaotic dynamics. Peaks in the propagator can be inter- 
preted as signatures of partial transport barriers within 
the chaotic sea. They arc visible as long as the parame- 
ter t of the velocity distribution is smaller than the time 
scale for crossing the barrier. As expected, for long times 
(t > 10^) only a narrow peak survives at a velocity which 
is in good agreement with the prediction of the sum rule 
(dashed line for t = 10^). 

Since the shape of the velocity distribution depends 
strongly on time, any definition of ballistic channels and 
the corresponding subdivision of the chaotic invariant set 
must be highly arbitrary. We single out the most promi- 
nent transporting islands which are visible in Fig. 4a close 
to the lower and the upper boundary of the chaotic sea. 
They have winding numbers W- ~ 30 and w+ ~ 20, re- 
spectively. In these regions we observe particularly sharp 
peaks in the velocity distribution for t < 10^ which are 
signatures of the corresponding Levy walks. Following 
Ref. [5] we continue by averaging the velocity distribu- 
tion over a region that excludes all such ballistic channels 
(-28 < V < +18 for the solid line in Fig. 4c; note that 
this t'-interval is defined with respect to the average ve- 
locity and therefore is not completely inside the chaotic 
layer in the Poincare section shown in Fig. 4a). The 
result represents the contribution from the bulk of the 
chaotic sea. It is definitely non-zero and in fact quite 
close to the asymptotic transport velocity (dashed line), 
irrespective of the time scale and the precise cutoff val- 
ues used. In other words, the average chaotic transport 
in this example is mainly due to the bulk region while the 
ballistic channels and their Levy walks contribute small 
corrections only. 

This shows that only the invariant sets, as featured in 
the sum rule Eq. (19), provide an appropriate concept 
for the description of the asymptotic directed transport. 



F. Biased ratchets 

Can Hamiltonian ratchets be used to transport parti- 
cles against an external force? As explained in Sec. II A, 
a constant force does not destroy the periodicity of the 
dynamics, and we can still resort to a unit cell to un- 
derstand the transport properties. The key question is, 
which invariant sets may survive in presence of an ad- 
ditional potential Vbias(a;) = ex. In Fig. 5a we compare 
two trajectories for c = 0.13 to the familiar phase-space 
portrait at c = (Fig. 2) . One of them was initialized on 
a large transporting island with winding number w = 1. 
Clearly, this island is still present although it is distorted 
and shifted in position. The winding number of the island 
is conserved since it is a topological quantity restricted to 
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rationals. Hence all trajectories inside the islands have 
asymptotic mean velocity iJ = 1 and we may conclude 
that Hamiltonian ratchets can transport uphill! This is 
confirmed by the full line in Fig. 5b which shows position 
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FIG. 5: (a) Stroboscopic Poincarc section for the model 
Eq. (5) as in Fig. 2. On top two trajectories of a system 
with the additional potential Vbias(a;) = cx with c = 0.13 are 
shown. One trajectory (big dots) was started at p = 10, i.e., 
in a phase-space region which is in the original system filled 
by non-contractiblc regular tori. In the presence of the bias 
such tori are absent and the trajectory keeps loosing momen- 
tum without bounds. In the extended system this trajectory 
is similar to a parabola (dashed line in (b)). The other tra- 
jectory (thick line in (a) and (b)) is part of a regular island 
with winding number w = 1, i.e., in the extended system 
this trajectory is transporting uphill without loosing momen- 
tum. The inset of (a) shows the shape of the regular island 
at different magnitudes of the bias potential. At c > 0.15 the 
island disappears. Also islands with negative or zero winding 
number do exist in the biased ratchet (not shown). 



vs. time for the same trajectory. 

The other trajectory was initialized in a phase-space 
region which for c = contains non-contractible KAM 
tori with positive winding numbers. We observe that for 
c = 0.13 the momentum of this trajectory is decreasing 
without bounds under the influence of the constant bias 
force, as naive expectation suggests. Only in a short 
time interval, when pt w 0, the driving potential has 
a relevant influence on this trajectory. For long times 
it behaves essentially like a free particle accelerated by 
the bias potential. Therefore x{t) for this trajectory is 
approximately parabolic (dashed line in Fig. 5b). 

From the presence of this single accelerated trajectory 
we can already conclude that no regular KAM tori sur- 
vive in the biased system (at least not in the phase-space 
region displayed in Fig. 5), since these would represent 
impenetrable barriers to transport in p-direction. Note 
that the KAM theorem does not apply to this situation: 
A constant force does not represent a smooth perturba- 
tion for the unit cell since the potential is not periodic. 
In fact there is a simple argument suggesting that an 
arbitrary small mean force destroys all non-contractible 
KAM tori: Assume that there is a KAM torus of the 
form p[^, t) periodic in ^ and r. Consider its average 
momentum at some given moment in time 

p{t)= r d^piCr). (25) 
Jo 

As we show by a straightforward calculation in Ap- 
pendix A the increment of p after one temporal period is 
given by 

p{t + 1) - p{t) = - [ dC dr V'{^, t) . (26) 
Jo 

Clearly, this increment must vanish for an invariant KAM 
torus. However, the r.h.s. of Eq. (26) is not zero for 
a biased system with a mean force. We conclude that 
no extended KAM tori survive and that therefore the 
chaotic sea is no compact invariant set anymore. Hence 
an arbitrarily small bias potential will destroy the chaotic 
ratchet transport in models like Eq. (5) while uphill 
transport can be realized by preparing initial conditions 
on regular islands of the phase space. 

G. A minimal model 

According to the previous sections, the decisive prop- 
erty of a Hamiltonian ratchet is an asymmetric mixed 
phase space. Based on this insight we can now construct 
minimal models for Hamiltonian ratchets which have this 
property and are otherwise as simple as possible. Prob- 
ably the simplest type of model with a mixed phase 
space are area-preserving maps generated from kicked 
one-dimensional Hamiltonians of the form 

H{x,p,t)=T{p) + V{x)J2S{t-n). (27) 

n 
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FIG. 6: (a) Poincarc section p vs ^ of a unit cell for the 
map given by Eq. (31). (b) Velocity distribution P{v) of 10'' 
trajectories started at random on the line p = 0, a: e [0, 1) in 
the chaotic sea of the system and iterated until 10^ and 10*, 
respectively. 



Integrating the equations of motion over one period of the 
driving we obtain an expUcit map expressing position x„ 
and momentum p„ immediately before the kick at t = n 
in terms of the values before the preceding kick 

Pn+l =Pn- V'{Xn) , Xn+1 = Xn + T' {pn+l) . (28) 

The most prominent example is the kicked rotor 

T{p) = Y ^(^) = ^ cos(27ra;) , (29) 

one of the best-studied paradigms of Hamiltonian chaos 
[20]. The phase space of this model is periodic with pe- 
riod I both in X and in p. Therefore one can define a 
compact unit cell with area Ax Ap = 1. 

The kicked rotor found an important experimental re- 
alization in the dynamics of cold atoms in pulsed laser 
fields [21, 22]. In this experimental setup the momen- 
tum instead of the position is the experimentally accessi- 
ble quantity and one is therefore interested in transport 
along the momentum direction. Apart from this purely 
formal difference, atom optics experiments promise to be 
ideal realizations of Hamiltonian ratchets. For this pur- 
pose one has to modify the phase space of the unit cell 



such that transporting islands arise and the symmetry 
x —X, p —>■ —p of the kicked rotor is destroyed. 

In fact transporting islands appear already in the stan- 
dard kicked rotor at kicking strengths K > 2iTm. They 
are referred to as "accelerator modes" [20] and leave 
traces in the dynamics which were also experimentally 
observed [22] . In the kicked rotor these accelerator modes 
always come in pairs transporting in opposite directions 
and therefore do not lead to transport in the chaotic sea. 
However, this symmetry can be destroyed, e.g., by ap- 
plying more than a single kick per period or by using 
asymmetric potentials in Eq. (27). It is not expected 
that the details of these manipulations will be of impor- 
tance for the resulting chaotic transport since, as we have 
shown in the previous sections, the latter is determined 
by the underlying phase-space structure only. 

In the remainder of this paper we therefore study an 
abstract model in the form of Eq. (28). The functions 
T{p) and V{x) are selected without reference to any par- 
ticular experimental setup and only guided by the desire 
to have a simple phase-space structure with a large trans- 
porting island. We choose 



V{x) = (ccmodl - 1/2)2/2 
T{p) = |p| -h3sin(27rp)/(47r2). 



(30) 



The resulting map 



Pn+l = Pn- {Xn mod 1) + 1/2 

Xn+i = x„ + Sgn(j;„+i) + 3cos(27rp„+i)/27r (31) 

is considered on a cylinder with transport along the ex- 
tended a;-axis while p = p+ 1 is here a cyclic variable that 
can be represented with p G [—1/2, +1/2). If the map is 
restricted to one unit cell x ^ ^ = x mod 1 we obtain the 
phase-space portrait shown in Fig. 6a. It shows one large 
regular island around the stable fixed point = 1/2, 
Po = —1/4 with winding number wo = —1. Due to the 
term \p\ in T(p) the phase space has no reflection symme- 
try around p = and also no other momentum-inverting 
symmetry such that there is no equivalent island trans- 
porting in positive direction. 

There are also no extended regular tori and the whole 
unit cell must be considered as the analogue of the com- 
pact stochastic layer in the continuously driven mod- 
els which we considered in the previous sections. Con- 
sequently the l.h.s. of the sum rule Eq. (21) vanishes, 
— WchVch + (— 1) Vreg- In Other words the total trans- 
port, averaged over the whole available phase space, van- 
ishes for this system which confirms that it is unbiased. A 
considerable simplification results from the fact that here 
the chaotic transport velocity can be computed from the 
relative phase-space volume of the single regular island 
Areg = 1 — ^ch alone. 



Veil — j4j-eg/ (1 j4j-eg) . 



(32) 



From the Poincare section Fig. 6a we find Aj-^g = 0.117± 
0.001, thus Vch = 0.133 ± 0.001. This is in very good 
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agreement with Vch = 0.1344 ± 0.0003 obtained directly 
from the spatial distribution of lO"' trajectories after 2 x 
10'' kicks. Fig. 6b shows the convergence of the chaotic 
velocity distribution to a delta function concentrated at 
this value, in accordance with Eq. (11). 

We would like to stress again that the directed chaotic 
transport in this ratchet model is a consequence of the 
phase-space structure and cannot be explained by the 
asymmetric kinetic-energy function alone. We have veri- 
fied this fact by repeating the analysis for a larger poten- 
tial 5V{x). Then the phase space is completely chaotic, 
yet despite the asymmetric function T{p) no transport is 
observed. 



III. QUANTUM RATCHETS 

We now turn to the investigation of quantized Hamil- 
tonian ratchets, i.e., driven ID Hamiltonian quantum 
systems which are classically periodic both in space and 
in time. We restrict attention to systems in which the 
phase-space volume of a unit cell is finite and phase space 
is composed of a chaotic sea with one or more embedded 
regular islands, as in the minimal ratchet model discussed 
above. This restriction leads to a finite Hilbert-space 
dimension which simplifies the numerical calculations. 
Moreover, we have seen that the dynamical processes rel- 
evant for transport are restricted anyway to the compact 
chaotic layer of the unit cell. We therefore expect mod- 
els with finite Hilbert-space dimension to capture also 
the essential features of quantized ratchet transport. 

A. Floquet operator and eigenstates 

For a system periodic in time, one can still construct a 
dynamical group with a single time-like parameter, which 
however now becomes discrete, measuring time in units 
of the period of the driving. It is generated by the unitary 
evolution operator over one period, 

U{t+l,t) = fcxp(^-^J^dtH{t)^, (33) 

where T effects time ordering. The computation of this 
Floquet operator is simplified considerably if H{t) is a 
kicked Hamiltonian as in Eq. (27). Then the time evolu- 
tion from time t = m — e immediately before the kick m 
to time t = m -\- 1 — £ immediately before the following 
kick can be expressed in terms of T{p) and V{x) as a 
product 

U = e-^np)/h^-ivix)/h (34) 

of two operators which are diagonal in the position or the 
momentum representation, respectively. The time evolu- 
tion of a state is obtained by successive multiplications 
by phase factors and fast Fourier transforms effecting a 
basis change. An additional simplification results if we 



consider p as a cyclic variable p = p + I, as is the case 
with the minimal ratchet model Eq. (31) to which our 
numerical results will be restricted. In this case the wave 
function is periodic in p with tp{p+l) = tp[p) and conse- 
quently the conjugate variable x is restricted to the dis- 
crete values Xn = nh. Here, h denotes the dimensionless 
ratio of Planck's constant to the phase-space area of the 
classical unit cell which we set to unity in Eq. (31). It is a 
well-known peculiarity of models with this property that 
the periodicity of the classical potential ^(a; -|- 1) = V{x) 
(or at least V'{x + 1) = V'{x) in the case of our minimal 
model) does not necessarily lead to a spatially periodic 
Floquet operator. The reason is that the potential is now 
restricted to discrete values Vn = V{xn) = V{nh) and 
periodicity is achieved only if there is an integer with 
Vn+ N — Vn which impHcs hN = M with another integer 
M. Hence h = M/N must be rational. In contrast, in 
periodic systems with infinite phase-space volume such 
as Eq. (5), the Floquet operator is spatially periodic ir- 
respective of the value of Planck's constant. In the fol- 
lowing sections we shall use values h = 1/N to ensure 
that the quantum system has the same spatial periodic- 
ity as the classical model. Only in the last Section III D 
we consider modifications of our results for irrational val- 
ues of h. They are to be interpreted as a spatial disorder 
that does not affect the classical phase-space structure 
but destroys the perfect periodicity of the corresponding 
quantum system. 

A double periodicity, both in space and time, requires 
to combine the corresponding representations of quantum 
mechanics appropriate for these symmetries, i.e., Bloch 
and Floquet theory, respectively. The eigenvalue equa- 
tion 

|.^„(t + 1)) = U\Mt)) = e-'""" \Mt)) (35) 

defines Floquet states |0c() and quasienergies £ [0, 1) 
[23] . For the systems considered here a is a discrete index 
l<a<N. 

For the discrete spatial translation group there is a 

continuous set of representations parameterized by the 
quasimomentum k £ [0, 1). In the simultaneous presence 
of temporal periodicity, the Bloch theorem now applies 
to Floquet states, 

cj,„,k{x + l,t)= e2'^"=<^a,fc(x, t) (36) 

so that both eigenstates and eigenphases carry a double 
index {a,k). The support of the Floquet band spectrum, 
in all cases considered here, consists of continuous lines 
in the two-dimensional {k, e)-space, cf. Fig. 7a. Since the 
spectrum is periodic with period 1 both in e and k, these 
variables are canonically conjugate to a pair of integers 
{nx,mt) which measure position and time in units of the 
spatial and temporal periods, respectively. We shall show 
in Section HI B that this allows to relate the band struc- 
ture to the time evolution of the spatial distribution by 
a double Fourier transformation. 

Wc have seen that the decisive property of classical 
Hamiltonian ratchets is the existence of invariant sets of 
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FIG. 7: (a) Quasienergy band spectrum of the minimal 
ratchet model Eq. (31) at = 32. Regular bands appear 
as approximately straight lines with negative slope, (b) The 
Husimi representation of the Floquet eigenstates correspond- 
ing to points on these lines are concentrated inside the regular 
island, (c) Most other eigonfunctions spread over the entire 
chaotic sea but avoid the regular island. The corresponding 
bands have strongly fluctuating slopes, (d) Distribution of 
band slopes (velocity expectation values) at = 128. The 
sharp peaJi sX v = —1 corresponds to the regular bands, the 
broader peak to the chaotic bands. The velocity of the clas- 
sically chaotic transport is marked by an arrow. 

phase space with different average velocities. Traces of 
the classically invariant sets are manifest in the quantum 
dynamics only if the quantum uncertainty allows to re- 
solve them, i.e., if h is much smaller than the relevant 
phase-space structures. From here on we shall restrict 
our attention to this semiclassical regime. Fig. 7a shows 
an example of a Floquet band spectrum for a Hamilto- 
nian ratchet, the minimal model (31). This system has 
two distinct invariant sets in phase space, the chaotic sea 
and one transporting island embedded in it. According 
to the semiclassical eigenfunction hypothesis [24, 25] one 
expects that in the semiclassical regime almost all eigen- 
functions condense on one of the invariant phase-space 
sets. Fig. 7b, c show the Husimi representations [26] of 
typical eigenstates. Indeed, one of them is concentrated 
inside the regular island while the other populates the 
chaotic sea, avoiding the island. Associated with these 
two types of eigenstates are two types of bands: regular 
bands appear in the spectrum as straight lines with slope 
dect/dfc ~ — 1, chaotic bands are fluctuating and have on 
average a positive slope. In the subsequent section, we 
are going to make this relation between bands and sub- 
sets of phase space more precise. We use it to establish 
a sum rule for transport in quantum ratchets analogous 



to the classical sum rule discussed in Section II D. 



B. Semiclassical transport in terms of Floquet 
bands 

1. Quantum sum rule 

The basic relation expressing the velocity of a Floquet 
state in terms of the quasienergy band to which it belongs 
is 

Voi,k = {{4>oc,k\'^\(l>oc,k)) = —TT' ■ i^'^) 

dk 

In the present case of a periodically driven system, the 
expectation value of the velocity operator i) = T\p) 
includes a time average over one period of the driving 
((...)) = Jq dt{...). The second member of Eq. (37) 
then follows from applying the Hellman-Feynman theo- 
rem, which was proven for time-periodic systems in [23] . 

A wave packet localized on the scale; of a single unit 
cell or narrower corresponds to a nearly homogeneous 
distribution in k. The corresponding mean velocity for a 
whole band a vanishes, 

K)fc = ^'dA:^ = 0, (38) 

as is implied by the periodicity of the bands. Averaging 
also over energy, i.e., summing over the bands, we find as 
velocity average over the total Hilbert space of the unit 
cell, 

Equation (39) can be considered the quantum- 
mechanical counterpart of the classical sum rule for trans- 
port, Eq. (21). Effectively, the quantum sum rule like 
the classical one refers to a finite subset of phase space. 
Here, the cutoff is introduced by the finite dimension of 
the basis used to span the Hilbert space of the unit cell 
in calculating the band spectrum. 

The crucial step for this quantum sum rule is the av- 
eraging along a given band a over the entire Brillouin 
zone, Eq. (38). In particular, this amounts to regarding 
all band crossings, however narrow, as avoided crossings. 
If k were considered a parameter with a fictitious time 
dependence, the quantum time evolution under a slow 
change of k would respect avoided crossings in exactly 
this manner. Therefore these bands are referred to as 
adiabatic bands [27]. 

It follows, conversely, that a finite mean velocity can 
be obtained if modified bands are constructed by con- 
necting band segments across all avoided crossings with 
a gap below some threshold. Such bands determine the 
time evolution under a fast change of k and accordingly 
are called diabatic [27]. They are not associated to a 
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fixed band index a and therefore need not be periodic 
in k. So for individual diabatic bands Eq. (38) does not 
apply, their mean velocity can be finite. We argue in the 
following that indeed it is diabatic bands, not adiabatic 
ones, which semiclassically correspond to invariant sets 
of classical phase space, and to which a relation between 
band structure and directed transport must refer. 

Fig. 7 provides numerical evidence to justify the assign- 
ment of invariant sets to diabatic bands. For example, 
the regular island with winding number — 1 is associated 
to straight-line segments in the spectrum, correspond- 
ing to a quantum velocity Va.k = —1 with very small 
fluctuations. In contrast, chaotic regions are represented 
by "wavy" band sections with strongly varying slope to 
which a precise velocity value cannot be assigned. In this 
sense, it is legitimate to talk of "regular" vs. "chaotic" 
diabatic bands. 

In the following we will reconsider the sum rule 
Eq. (39) using diabatic bands and express the different 
contributions in terms of the invariant sets of the classical 
phase space. First we note that replacing in Eq. (39) adi- 
abatic by diabatic bands amounts to interchanging band 
indices at avoided crossings, thus it results at most in a 
permutation of terms within the sum but does not af- 
fect the sum rule as a whole. We can therefore group 
diabatic-band terms in Eq. (39) according to the classi- 
cal invariant set they pertain to. 



= 
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(40) 

In the semiclassical limit the respective numbers of terms 
in the sums are given by the relative fraction of phase 
space occupied by the corresponding invariant sets, i.e., 
-^ch = /ch N for the chaotic bands and = frN for the 
various embedded regular islands r. N = is here the 
total number of bands, i.e., the Hilbert-space dimension 
per unit cell. Assuming that the classical phase space 
contains only a single chaotic component we can char- 
acterize the associated diabatic bands by a mean slope 
(vch) and have J2aech{'^^a.k/<i-k)k = A^ch(wch). 

For the regular bands, the double periodicity of the 
(fc, e)-space allows to define winding numbers in the same 
way as we did in Section II B for the topology of regular 
islands in the conjugate (x, t)-spacc. For the same rea- 
son as the classical winding numbers these topological 
quantum numbers have to be rational, i.e., w^™ = n/m 
if the band closes upon itself after n revolutions in e- and 
m revolutions in fc-direction. As the regular states are 
localized on the invariant tori inside the island, their ve- 
locity expectation (band slope) in the semiclassical limit 
approaches the regular transport velocity. This leads to 
the conclusion 



der 
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(41) 



Avoided crossings modify the band slopes in a range 
which is negligible in the semiclassical limit (see 



Sec. IIIC), while the winding numbers as topological 
quantities are not affected at all. In other words, the 
winding number wJJ™ of a diabatic band r pertaining to 
a classical regular island r is identical to the classical 
winding number w'j} of that island, Eq. (8). We have 
now 



= Nf,i,{v,i,)+Nj2frvf- 



(42) 



Note that fch, fri and v^' are all classical quantities. 
Consequently, also the quantum transport velocity {vch) 
must coincide with its classical counterpart 



(43) 



This is the main result of the quantum mechanical sum 
rule. We stress again that it pertains to the semiclassical 
regime since otherwise the notion of diabatic bands is not 
applicable. 

Fig. 7d confirms Eq. (43) qualitatively. It shows the 
distribution of quantum velocities (band slopes) for our 
minimal ratchet model. We observe two well separated 
peaks, one for the regular bands at = —1 and one 
at v'^l^ for the chaotic bands. The region separating the 
two peaks corresponds to the band slopes in the vicinity 
of avoided crossings between regular and chaotic bands. 
The weight of the distribution in this intermediate region 
decreases with h and vanishes in the semiclassical limit 



2. Form factor 

Our analysis based on winding numbers can be ap- 
plied to predict the mean quantum transport velocity in 
the semiclassical regime from the classical value. The 
band spectra, however, contain more detailed informa- 
tion about quantmn transport, encoded in the spectral 
two-point correlation functions. A double Fourier trans- 
form e ^ rut, k ^ Ux and subsequent squaring of the 
spectral density translates two-point correlations in the 
bands into the entire time evolution of the spatial distri- 
bution on the scale of the temporal and spatial periods, 
respectively. 

As a suitable quantity to establish this relation, we re- 
cur to the generalized form factor introduced and studied 
in [28] for completely chaotic systems. We define it as 



K{nx,mt) = —{\u{nx,mt)\ 



(44) 



with 



u{nx,Tnt) 
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N 

= ^Ua{nx,mt). (45) 
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denotes the Hilbert-space dimension per unit cell, that 
is the phase-space area of a unit cell in units of Planck's 
constant h. Uk is the N x N Floquet operator (33) eval- 
uated at Bloch number k. The integers n^, rrit are the 
discrete variables canonically conjugate to k and e, re- 
spectively, that is, the unit-cell number relative to the 
starting point, and time in units of the period of the driv- 
ing. The average (...) in Eq. (44) is essential in order to 
remove the otherwise dominant fluctuations around the 
mean value. It can be taken over a narrow time range or 
over an ensemble of quantum systems corresponding to 
approximately the same classical system. 

As we will now show, the form factor is related, on 
the one hand, to the classical dynamics of a distribution 
which initially covers homogeneously the phase space of a 
single unit cell. On the other hand, it contains the quan- 
tum velocity distribution as a limiting case. Therefore it 
is an appropriate starting point for a semiclassical theory 
of ratchet transport. 

We assume to be sufficiently close to the semiclassi- 
cal limit » 1 such that we can consider the band 
spectrum in the diabatic approximation. Moreover, in 
the semiclassical limit it is justified to neglect correla- 
tions between diabatic bands pertaining to different in- 
variant sets (regular or chaotic) unless they are related 
by symmetries. This allows to write the form factor as an 
incoherent sum of the respective contributions, because 
the averaging in Eq. (44) suppresses uncorrelated cross 
terms. We obtain 

K{n^, mt) = ^ Krin^, mt) + Kchin^, rrit) , (46) 

r 

the sum running over all regular invariant sets (islands 
and island chains). 

In Appendix B we obtain the semiclassical expression 

Kr{nx, mt) = fr lJir5n^n^-Vrmt ■ (47) 

for the form factor of a chain of regular islands with wind- 
ing number Wr = Vr/ (Jr- It seems that the form factor is 
enhanced by a factor fij. for an island chain as compared 
to a single island of equal total size, but this is not the 
case. In Eq. (47), S^^n^-Vriyn = 1 holds only at the unit 
cell Tlx = {vr/ ^J'r)mt = VrTTit to which a classical trajec- 
tory, started in the regular island at n^; = 0, has traveled 
in time m*. In particular, as Ux is an integer, mi must 
be an integer multiple of /i^. That is, Kr{nx,rni) is fi- 
nite only every fi^-ih. period of the driving, such that the 
average contribution to the form factor is independent of 
the period Hr of the island chain. 

For the chaotic contribution to the form factor we can 
resort to a semiclassical theory which has been developed 
for completely chaotic systems in [28, 29]. In order to 
apply it to a system with a mixed classical phase space 
we assume the validity of the ergodic sum rule [30] for the 
chaotic component. Then the result of [28, 29] remains 
essentially unchanged, and the form factor is given in 
terms of the classical velocity distribution of the chaotic 



component as 

-ft^ch(nx,mt) = — Pch ( — ,TOt ) {mt<m]i). 

(48) 

To be precise, Pch(w,t) entering this equation is the 
chaotic classical propagator for a uniform distribution 
inside the chaotic introduced in Section HE. Its 

definition is Eq. (9) with po = Xch- Since Eq. (48) is based 
on the diagonal approximation [31], i.e., correlations be- 
tween different classical orbits have been neglected, it is 
valid only for short times and breaks down beyond the 
Heisenberg time mn ~ -/Vch ~ /ch-^ of the chaotic com- 
ponent. 



3. Quantum velocity distribution 

A complementary approximation to the form factor for 
long times can be achieved following again Refs. [28, 29] . 
The chaotic bands fluctuate as a function of k with an 
amplitude approximately given by the spacing Ae « 
between them. For times beyond the Heisenberg time, 
these fluctuations give rise to phase oscillations in the 
integrand of Eq. (45) which exceed 27r. Therefore we can 
perform the fc-integration in stationary-phase approxima- 
tion and obtain 

Ua{nx,mt) = V^/RU^ti 

X exp(27ri[A:snx - Ca.fesmt]) , (49) 

i.e., only those points k = ks contribute to the integral 
where the derivative of the phase of the integrand van- 
ishes, = rix — j.^mt. These are isolated points in the 
spectrum which can be assumed to vary independently 
upon averaging in Eq. (44). Therefore we can neglect 
all cross terms when squaring the sum of contributions 
from different points of stationary phase and obtain for 
the form factor 

Now that we are rid of all phase factors it is very instruc- 
tive to rewrite the result again as an integral over the 
Bloch number k 

i^ch(n„mO = ^— ^ I'dkdfe'^^k-—] ■ (51) 
mnmt^jQ \ ' rritj 

This equation has two important consequences. First we 

note that up to normalization the form factor beyond the 
Heisenberg time is nothing but the distribution of band 
slopes alias quantum velocities 

K{nx,mt) ~ Pqua.nt{v)\v=n^/mt ("^t > "^h) , (52) 
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which is shown for the minimal model in Fig. 7a. As in 
the classical case, this velocity distribution is the natu- 
ral quantity to describe a system with directed ballistic 
quantum transport and the form factor can be considered 
a useful generalization of it. 

Secondly, Eq. (51) implies that the form factor at any 
time mt beyond the Heiscnberg time mn can be ex- 
pressed via scaling by the form factor right at the Heisen- 
berg time 



mt V "^t 

"ch — , toh [mt > mn) (5c 

mt \mt J 



3) 



In the second line we have used the semiclassical approx- 
imation Eq. (48) for mt = mu- It is valid only up to the 
Heisenberg time, but according to Eq. (51) it determines 
the form factor also beyond. Of course, the validity of 
Eq. (53) depends on applying both, the short-time and 
the long-time approximations for the form factor right at 
the Heisenberg time where they arc on the verge of break- 
ing down. This interpolation procedure has been corrob- 
orated by comparison to results from numerics and from 
supersymmetry in [28, 29]. We expect that it applies in 
the present case of a transporting chaotic component as 
well. 

The two consequences of Eq. (51) combine to the con- 
clusion that the distribution of quantum velocities in the 
chaotic component of the band spectrum is equal to the 
distribution of time-averaged classical velocities for an 
ensemble of particles filling the chaotic component of 
phase space homogeneously. Information on the quan- 
tum system enters into this classical distribution only 
via the point in time at which this velocity distribution 
is evaluated it must be chosen as the Heisenberg time 
iVch of the chaotic component. Before writing down this 
result we note that the restriction to the chaotic com- 
ponent is actually not necessary, since for the embedded 
regular islands the same result applies trivially because 
of Eq. (41). Hence we have 



P. 



quant 



(v) = Pclasa{v,mu) , 



(54) 



for a stochastic layer including one chaotic component 
and one or more embedded regular islands. Eq. (54) is a 
nontrivial result because it establishes quantum-classical 
correspondence for the velocity distributions and thus 
for asymptotic long-time transport properties. We stress 
again that this result was derived semiclassically within 
the diagonal approximation. It would be very interesting 
to explore possible corrections due to neglected interfer- 
ences between classical periodic orbits (akin the weak- 
localization correction in the standard form factor [31]), 
but at present the methods to deal with such corrections 
[32] are not sufficiently developed to treat the type of 
system we are dealing with here. 



C. Long-time quantum transport and dynamical 
tunneling 

1. Transport of wave packets 

So far we have considered transport only in terms of 
stationary quantities like eigenstates and band spectrum. 
Using the obtained results we can now describe the trans- 
port of arbitrary wave packets. The asymptotic quantum 
transport velocity of a wave packet is an average over all 
band slopes, weighting each Floquet state by its over- 
lap with the initial state. To see this we write the wave 
packet as a superposition of Floquet states 
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X,t) = y2 <^klpa,k{t)'Poc,k{x) 
a Jo 

= V/ dfc^„,fee-2---'=*</,„,fc(a;), (55) 
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calculate the expectation value {x{t)) of position as a 
function of time (see App. C), and obtain 



r' + oc 



o{t)) = I dxx\tlj{x,t)f 

= Voot + 0{t) , 



with 
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(56) 



(57) 



Consider now a wave packet localized initially within 
a single unit cell and, inside this unit cell, on one of 
the invariant sets of the classical dynamics. Then the 
weights IV'q,*;!^ are approximately homogeneous in k but 
concentrated on the diabatic bands corresponding to the 
supporting invariant set. This is illustrated in Fig. 8. 
Consequently, the asymptotic velocity is an average over 
the corresponding band slopes. For example, for a wave 
packet started inside the chaotic sea we expect a value 
close to the classical chaotic transport velocity because 
this is the average slope of the chaotic bands, see Eq. (43). 
We confirm this semiclassical result in Fig. 9, where the 
average position of two chaotic wave packets is shown 
over a large time interval and for two different values of 
N = h~^. In agreement with Eq. (56), we observe a lin- 
ear dependence on time with very small fluctuations, i.e., 
asymptotically there is indeed directed ballistic quantum 
transport. The precise value of the velocity depends on 
the initial conditions but these fluctuations decrease with 
h and the average approaches the classical transport ve- 
locity. Typically the quantum velocity for a semiclassical 
chaotic wave packet is slightly above the classical value. 
This is a consequence of the hierarchical phase-space re- 
gions around the embedded islands which communicate 
with the main chaotic sea only via leaky cantori. Depend- 
ing on h, quantum transitions across some of these can- 
tori are possible only by tunneling, i.e., they are almost 
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FIG. 8: Quasienergy band spectrum of the minimal ratchet 
model at = 32. The linewidth encodes the overlap 

Ki?^a,fc|V')r of the corresponding Floquet state 10a, fc) with an 
initial wave packet In (a) this wave packet is a coherent 
state located in the chaotic part of the phase space of a single 
unit cell, in (b) it is concentrated on a torus inside the major 
regular island (cf Fig. 6a). 



blocked. Therefore the part of the chaotic component en- 
closed by these cantori effectively belongs to the regular 
island [33] and, according to the sum rule (21), this en- 
hanced island size is compensated by a correspondingly 
larger chaotic transport velocity. 



2. Dynamical tunneling 

On first sight it may surprise that the division of clas- 
sical phase space into invariant sets can influence the 
long-time quantum dynamics. After all, classically im- 
penetrable barriers can be crossed in quantum dynamics 
by tunneling. Tunneling is known best for the case of en- 
ergetic barriers, e.g., in a double- well potential. Dynam- 
ical tunneling is the generalization of this phenomenon 
to barriers in phase space [34] and was recently demon- 
strated experimentally [35, 36]. If in quantum dynamics 
no strict barriers exist, the wave packet should explore 
the entire accessible phase space for sufficiently long time 
and consequently directed transport should vanish, at 
least on average. Wc have seen in the previous section 
that this is not the case. So what is the role of tunneling 
in Hamiltonian ratchets? 

To answer this question we consider a wave packet 
which is initially prepared inside the regular island within 
the unit cell = 0. Classically, such an initial distribu- 
tion is simply transported along the chain of regular is- 
lands with a velocity corresponding to the winding num- 
ber Wr, i.e., Pr{x + wvt, t) = Pr{x, t = 0). This property 
is conserved in the quantum dynamics if we neglect the 
narrow avoided crossings in the band spectrum which ac- 
count for the difference between adiabatic and diabatic 
bands. Let us demonstrate this for the regular island in 
our minimal model which has winding number Wr = — 1- 
The diabatic regular bands are straight lines with slope 



I.e., 



fc = Er + Wrk . 



(58) 



As illustrated in Fig. 8b, a localized initial wave packet 
can be constructed from such a band by a uniform su- 
perposition of all states 
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FIG. 9: Position vs time for wave packets initialized as co- 
herent states inside the chaotic part of phase space of the 
minimal ratchet model Eq. (31). Two different values of h 
and two different initial conditions are used. The dotted line 
shows the classical chaotic transport velocity. 
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We restrict attention to times which are a multiple of the 
period Hr of the central orbit inside the island. Then Wr t 
is an integer which indicates one particular unit cell. We 
measure x relative to that unit cell and find for the wave 
packet 



'^{X + Wrt, t) 
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dk exp(— 27rier,fet) 4'r,k{x + Wr t) 



dk exp(27ri(A:Wr — ^r,k)t) 4'r,k{x) 



exp(-27ri e,.,oi) *(.t,0) , 



(60) 



This shows that the wave packet is indeed transported 
like the corresponding classical distribution. It has the 

asymptotic velocity Wj. and docs not show any spreading, 
i.e., there is no signature of dynamical tunneling within 
the approximation of diabatic bands. 

We conclude hat tunneling out of an island in classical 
phase space is encoded in the avoided crossings between 
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the regular and the chaotic bands. These avoided cross- 
ings show up in the regular bands as deviations from the 
straight line 6^,0 + ijJvk- Close to an avoided crossing the 
regular bands are bent towards the chaotic bands, i.e., 
the actual slope is fc-dependent and slightly smaller than 
w,. . Using this qualitative information about the shape of 
the regular bands we can make a prediction for the shape 
of the wave packet at very large times t —>■ oo. In this 
regime the wave packet can be calculated from Eq. (55) 
in stationary phase approximation. We find 
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/ dk exp(27ri(fcX — Cr^kt)) (l>r,k{x) 
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X exp(27ri [fce^ ^. - er,k]t) (pr,k{x) . 



(61) 



We have again decomposed position into a large inte- 
ger X denoting the unit cell and the remaining fraction 
< a; < 1. (j>r.k{x) is considered a slowly varying pre- 
factor of the rapidly oscillating phase. The points of sta- 
tionary phase in Eq. (61) select the Bloch states whose 
superposition yields the wave packet at time t and po- 
sition X. It is no surprise that these are exactly the 
points for which the slope of the band corresponds to 
the velocity X/t. Due to avoided crossings, the actual 
slope of the regular bands is smaller than Wr- Hence for 
the transition to the unit cell X = w^t where all clas- 
sical probability is concentrated, no points of stationary 
phase with real k exist: To leading order this process is 
forbidden in quantum mechanics! There might be com- 
plex solutions of the equation e'^ ^ = Wr, but then the 
exponent in Eq. (61) has a real part and the contribution 
will be exponentially small in t, which is indeed observed 
in Fig. 10(b). The main part of the wave packet is con- 
centrated not in the "classical" unit cell but rather at 
positions for which real points of stationary phase exist 
in Eq. (61). These correspond to velocities distributed 
narrowly around a value slightly below the classical ve- 
locity. Due to this dispersion in the velocities, induced 
by avoided crossings, the wave packet will spread ballis- 
tically in time and will be peaked behind the classically 
expected position (Fig. 10(a)). 

For a wave packet initially prepared in the chaotic 
part of a unit cell the influence of tunneling is much less 
pronounced (not shown): Although the narrow avoided 
crossings with regular bands do modify the chaotic bands 
as well, the existence of points of stationary phase in an 
expansion similar to Eq. (61) is unaffected: due to the 
wide avoided crossings betwecin themselves, the chaotic 
states have a large variation in their velocities around the 
classical value anyway. 

We have thus identified the role of tunneling in Hamil- 
tonian ratchets. It leads to avoided crossings between 
regular and chaotic states (or between regular states 
with different winding numbers) . In the dynamics of ini- 
tially localized wave packets tunneling shows up mainly 
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FIG. 10: (a) Black line: Wave packet prepared in the regular 
island of the unit cell a; = and propagated to time t = 1000 
in the minimal ratchet model Eq. (31) with = 16. 

The classical probability would be restricted to the unit cell 
X — —1000, while the quantum wave packet has tunneled out 
of this "classical" unit cell and starts spreading. However, 
there is a large peak lacking slightly behind the classically ex- 
pected position. Gray line: Same for irrational = 16 + a. 
In this case the Floquet operator has no spatial periodicity. 
The part of the wave packet outside the classical unit cell 
localizes and develops an asymmetric envelope with approx- 
imately exponential tails. Inset: The probability to remain 
inside the classically expected unit cell x = —1000 is the same 
for rational and irrational h. (b) Due to dynamical tunneling 
the quantum probability in this "classical" unit cell decays ex- 
ponentially as a function of time. With respect to this decay 
the periodic model with h^^ — 16 is almost indistinguishable 
from the aperiodic model with irrational = 16 -|- cr. 



in the evolution of regular states, which slightly lag be- 
hind the position expected from classical considerations. 
We stress again that tunneling is not able to hinder di- 
rected ballistic transport of such wave packets even for 
infinite time. 

An interesting and important special case are systems 
with a symmetry-related pair of counter-moving regular 
islands like the kicked rotor in presence of accelerator 
modes. Dynamical tunneling between such island pairs 
was demonstrated experimentally [35, 36]. It is crucial 
to understand the difference between our argumentation 
above and this situation. First we note that a pair of 
symmetry-related islands is not analogous to a symmet- 
ric double-well potential. In the latter case all eigenstates 
are superpositions of left and right. Below the barrier 
top, their eigenenergies form quasidegenerate doublets 
and thus contribute to tunneling. In the case of counter- 
moving islands this applies only to the vicinity of avoided 
crossings between the corresponding bands where indeed 
they form a doublet. Away from these isolated and semi- 
classically small regions in fc-space the bands are approxi- 
mately straight lines but with opposite slopes, i.e., there is 
no systematic degeneracy. In this paper we consider wave 
packets initially localized inside one unit cell. In A;-space 
such a wave packet is extended. Therefore its weight in 
the vicinity of avoided crossings, where it contributes to 
tunneling, is negligible. By contrast, in the experiments 
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mentioned above the wave packets extend initially over 
many unit cells. Therefore, in fc-space they may well be 
concentrated right at avoided crossings. Then, and only 
then, dynamical tunneling is the expected consequence. 

D. Quantum transport in the presence of disorder 
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FIG. 11: (left) Averaged velocity expectation values of 
wave packets initialized in the classically chaotic region at 

(a;o,po) = (0, 0.25 ± 0.05). Beyond tioc, the velocity oscillates 
around zero, (right) Time tmax at which the averaged velocity 
expectation value falls short of a given limit under variation 
of h. For numerical reason, we chose two different values for 
■yiim. 

In this last section we will describe some modifications 
of the quantum transport in a situation, when the exact 
quantum periodicity is destroyed by weak static disorder. 
As explained above in Section III A this can be realized 
easily within our minimal ratchet model by choosing an 
irrational value of h. In this case the Bloch theorem does 
not apply anymore and, on a large scale, we expect dy- 
namical localization of wave packets and eigenstates. The 
properties of the eigenstates and in particular the failure 
of the semiclassical eigenfunction hypothesis in this case 
have been studied in [15]. We will here concentrate on 
the evolution of wave packets in the presence of disor- 
der. In Fig. 10 (gray line) we display the shape of a wave 
packet which was initialized in the regular island of unit 
cell X = a.t time t = 1000. Initially the wave packet 
follows the classical evolution, i.e., it is transported at 
velocity v = ~1 and looses probability due to tunneling. 
The process of timneling out of the island is essentially 
the same as in the case of a periodic system with ra- 
tional h. This is demonstrated by Fig. 10b and also by 
the inset of Fig. 10a, where one can see that the proba- 
bility remaining inside the classical unit cell is the same 
for both systems. However, the fate of the probability 
which has tunneled out of the island is entirely different 
from the periodic case. We see in Fig. 10a that the wave 
packet develops exponential tails which are characteristic 
of localization. Unlike the periodic case, the max;imum 
of the wave packet is not close to the classical expec- 
tation but rather close to the origin, i.e., the disorder 



prevents quantum transport despite the underlying clas- 
sical ratchet mechanism. The latter is manifest, however, 
in the asymmetric shape of the wave packet which has a 
much longer tail in the direction of classical transport. 

Similarly, disorder does also affect wave packets which 
arc initialized in the chaotic sea. Fig. 11a shows the ve- 
locity expectation value for such a wave packet at two dif- 
ferent values of the effective Planck's constant h. There 
is an initial period when (v) ^ Vch, but then the veloc- 
ity drops to zero because the wave packet tunnels into 
the island and finally occupies the whole available phase 
space. The time for this process is expected to scale as 
t ~ e'^/'' [37]. As Fig. lib shows, this is also the time 
scale for which the quantum ratchet shows transport in 
the presence of disorder. This maximum ratchet opera- 
tion time tniax can be defined as the time at which the 
velocity of a wave packet falls below a certain threshold. 
In Fig. lib log tmax is seen to be approximately propor- 
tional to h~^. Hence, in the deep semiclassical regime 
the quantum ratchet can work over an exponentially long 
time even in the presence of static disorder. 



IV. DISCUSSION 

The study of ratchets has largely been motivated by 
the interest in the physical principles of intracellular 
transport: Motor molecules, driven by chemical energy, 
are moving along chain molecules whose length is of the 
order of the cell size, and which consist of millions of units 
concatenated in a highly ordered manner, resembling the 
crystal order encountered in inorganic solids. It is there- 
fore natural to model them as one-dimensional, infinitely 
extended potentials with exact spatial translation invari- 
ance, but with reflection symmetry manifestly broken to 
define a preferred direction of transport. 

While the breaking of mirror symmetry is crucial to 
obtain directed transport, the role of translation invari- 
ance appears circumstantial, at most of heuristic impor- 
tance for the theoretical description. Translation invari- 
ance has been indispensable, however, in order to achieve 
first analytical and numerical results on directed trans- 
port in ratchets. In the present context of Hamiltonian 
systems, it allowed us to show that directed transport 
comes about by counter propagating phase-space flows 
within regular and chaotic components of systems with a 
mixed phase space. Moreover, quantum ratchets are ob- 
tained by quantizing Hamiltonian ratchets in the frame- 
work of Bloch theory; they exhibit transport at similar 
rates as their classical counterparts, at least in the semi- 
classical regime. 

Real systems showing directed transport, biological or 
physical, though, break translation invariance in various 
ways and to various degrees, the only exception being 
systems where the spatial coordinate is cyclic, like in bi- 
ological "rotation motors" or in pumping devices in a 
closed configuration [38]. In the following we discuss a 
number of typical deviations from spatial periodicity and 
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their consequences for transport. Since quantum systems 
are far more sensitive to the presence or absence of sym- 
metries than classical ones, the question concerning im- 
perfections of translation invariance becomes even more 
crucial on the quantum level. 

Experimental realizations of Hamiltonian ratchets, as 
in optical lattices or in solid-state devices, always show a 
certain amount of disorder, in the form of small stochas- 
tic differences between the unit cells. Classically, smooth 
spatial disorder, if it is not too strong, will not completely 
disrupt the phase-space structures underlying transport 
on short time scales, thus it has only a minor effect on 
transport [39]. For long times, however, we expect that 
transport is destroyed. In extended quantum systems ar- 
bitrarily weak randomness in the potential immediately 
leads to localization. As we show in Section HID, even a 
type of disorder that is invisible in the classical dynamics 
entails a breakdown of quantum transport on a timescale 
proportional to the localization length. It should be 
kept in mind, however, that localization as a quantum 
coherence effect is counteracted, in turn, by incoherent 
processes caused by the unavoidable coupling to ambi- 
ent degrees of freedom, or similarly by a "noisy" driving 
that breaks temporal periodicity. While it is well known 
that in this way, incoherence partially restores diffusive 
transport in systems with dynamical localization [40] , its 
effects on directed transport remain to be explored. 

The presence of a spatially homogeneous force breaks 
translational invariance in a more controlled yet radical 
manner. Rather than forming an unavoidable nuisance, 
it may be imposed intentionally to extract work from a 
ratchet. Moreover, it allows to define a stall force as 
the external bias just sufficient to bring transport to a 
standstill [41], and to ascribe an efficiency to ratchets. 
In contrast to disorder, a finite mean potential gradient 
forms a perturbation of unbounded amplitiide, and thus 
radically changes the structure of the classical ratchet 
phase space. Still, as explained in Section II F, directed 
regular transport reacts smoothly on an external bias, 
i.e., it requires a gradient of the order of those present 
in the original periodic potential to be completely sup- 
pressed. On the quantum level, additional complications 
arise in that eigenstates become metastable and eigenen- 
ergies correspondingly complex. This situation can be 
handled in a framework similar to scattering theory [42] . 
Its application to ratchets is under way. 

Finally, in most physical setups, transport takes place 
between two "terminals" , typically modelled as electron 
reservoirs. This amounts to confining the ratchet proper 
to a finite section of space- yet another elementary way 
to break translational symmetry. Taking it into account 



would allow to make contact with a different, but closely 
related paradigm of directed transport: Pumps are de- 
vices that channel a well-defined amount of charge, mass, 
etc., per cycle of an applied force from one terminal to 
the other [43, 44]. Obviously, pumps can be considered 
as ratchets reduced to a finite number of unit cells, or 
conversely, ratchets could be constructed by concatenat- 
ing an infinite number of pumps or equivalently, by clos- 
ing the pumping circuit. The only difference lies in the 
kind of model usually studied in these respective con- 
texts, namely fast drivings resulting in a chaotic dynam- 
ics in one case, slowly driven potential wells that resemble 
peristaltic pumps in the other [43]. But this is an artifi- 
cial distinction: It has been shown recently that driven 
chaotic scattering systems, employed as pumps, also gen- 
erate directed transport if all relevant binary symmetries 
are broken [45]. 

In order to study ratchets as realistic devices clamped 
between reservoirs at given temperatures and chemical 
potentials, however, another crucial building block is 
missing, a quantum statistical theory of transport under 
strong time-dependent driving far from equilibrium. For 
first approaches to this problem from the points of view 
of quantum scattering and quantum transport theory, 
see Refs. [46] and [38], respectively. 
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APPENDIX A: CHANGE OF MEAN 
MOMENTUM OF A KAM TORUS 

In this Appendix we consider a non-contractible KAM 
torus that can be specified by the functional dependence 
of the momentum on position and time p{^, t). Note that 
the existence of such a function is an assumption which 
simplifies our reasoning. 

Wc consider the average of the function p{(_, t) along 
the torus and replace the integral representation of this 
quantity, Eq. (25), by a Riemann sum over N ^ oo dis- 
crete points = n/N, Pn = p{^n, r) 



N 



(Al) 



Similarly we introduce a discrete time increment 5t and 
find that the A*" phase space points {^mPmT) in Eq. (Al) 
evolve to {S,mPn, t + 6t) with 

in=^n+ PnSr Pn = Pn - V'{^n,r)ST . (A2) 

Now wc use these new points to discretizc the integral 
representing p(t + St). In this way we obtain an expres- 
sion for the time derivative of p, which we evaluate to 
leading order in St and and then transform back to 
an integral. We obtain 



dr 



p{i 



1 ^ - 

n— 1 

= S^J2{lN + ^r.+l-Pn)ST 
n=l 

x[pn-V'{^n,r)ST -j^Pn] 

~ J^^'(^n,T) + {Pn+l-Pn)Pn] 

n=l 

^ j;fT.{-^'(^-,r)+PnP'{^n,r)} 

n=l 

» / d«{-V"({,T)+P«,T)^pK,T)} 

= /<«{-n«,x) + i|p"«,.)} 



/ 

Jo 



d^V'{^,T) 



(A3) 
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For the last line wc have used the periodicity of the func- 
tion r) with respect to ^. By integration with respect 
to T we find Eq. (26) which was the purpose of this Ap- 
pendix. 



APPENDIX B: GENERALIZED FORM FACTOR 
FOR AN ISLAND CHAIN 

In this appendix we derive Eq. (47). We consider 
the contribution to the form factor from one particu- 
lar chain of regular islands r. If the winding number 
is wf = Vr/ l^r then inside a unit cell this island chain 
consists of islands which are traversed in sequence. 
In the somiclassical limit, wc associate (diabatic) bands 
with index a to the island chain. These bands consist 
of straight line segments with the slope = wf, cf. 
Eq. (41). The segments are connected such that the dia- 
batic band as a whole is periodic in e and fc, with periods 
Vr and /i,., respectively. It is easy to see that for a given 
value of k there are jir equidistant segments (values of 
the quasienergy) pertaining to the same diabatic band 
a. Semiclassically, the number of states associated to 
the island chain for given k is approximately frN where 
fr is the fraction of phase space occupied by the island 
chain as a whole and N = hr^ is the total number of 
bands. It follows that the number of complete diabatic 
bands associated with the island is frN/ jir- 

To integrate over a diabatic band consisting of many 
straight segments it is convenient to consider instead an 
extended Brillouin zone in which the band corresponds 
to a single straight line 

er,a,fe = ( er,a,o + —k \ mod 1 , fc G [0, iir) . (Bl) 

In this way we can perform the /c-integration in Eq. (45) 
and find 

Ua{n^,mt) = r''dA;e['"'('=»--(''-.=.o+7^'=)'»*)] 

^0 

Jo 



-27rier 



fir e-2'^''-°'°'"'(5, 



(B2) 



For the contribution of the island chain r to the form 
factor we have now 



(B3) 



a=l 



i.e., wc have to perform a sum over quasicnergies at fixed 
Bloch number A: = which can be done in the same way 
as for the spectrum of eigenenergies pertaining to regular 
states of an autonomous system [47, 48] . We assume the 
dynamics within the island to deviate sufHciently from 



harmonic vibrations around its central orbit. Then the 
spectrum of quasicnergies er,Q,o will not be equidistant 
and the phases in Eq. (B3) from different a can be as- 
sumed uncorrelated in the semiclassical limit. This allows 
to replace | J2a • • • P number of terms in the sum, 

which finally yields Eq. (47). 



APPENDIX C: WAVE PACKET TRANSPORT 

We compute the average position of a wave packet 

ijj{x,t) for long time t ^ 1. First we write the wave 
packet as a superposition of Floquet eigenstates (j)a,k {x) 
with quasienergy €a,k 

1p{x,t) = dk'llJa,k{t)'Pa,k{x) 

= y2 I cifcVa,fce-2---'=*<^„,fc(x), (CI) 

a Jo 

where 

/ + 00 
dxCfe(a^)^(a^,i = 0). (C2) 
-oo 

The integral representing the expectation value of x for 
the wave packet (CI) can be split into two contributions 
^, X corresponding to length scales within a unit cell and 
over many unit cells, respectively 

/+cx> 
dxx\%l;{x,t)\^ 
-oo 

r.1 +00 



= dx {x-\-n) \t{j{x -\- n^t)\^ 
Jo 



n— — oo 



(C3) 



Naturally, the contribution from the dynamics inside the 
unit cells is bounded from above by the size of the unit 
cell 

C(t) = / dxx J2 \i^{x + n,t)f 

Jo n=-oo 

< dx \^{x + n,t)f 

Jo „=-oo 

= 1 (C4) 

(the last equality expresses the normalization of the wave 
packet). Therefore ^ is irrelevant for directed ballistic 
transport. 

Evaluating the term that describes the wave packet 
on large scales, we use the Bloch theorem to switch from 

position representation to the conjugate variable k, where 
a spatial shift corresponds to differentiation. We have 



nil){x + n,t) 



dkipa,k{t) (l>cc,k{^ + n) 
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= J2 [ dfcVa,fc(i)</'a,fe(a^)rie'"'*'" 

a -^0 



1 ^2Trikn 



^ZTTlKn ^ 

27ri dfc 



The last line follows from partial integration and the pe- 
riodicity in k of e^'^''^"i/'a,fe(t) 4'a,k{x)- Inserting this into 



X{t) = dx ^ n\ip{x + n,t)\ 
Jo 



n— — oo 



and decomposing also the complex conjugate ip*{x + n, t) 
into Floquet states we find 



X{t) = - dx J2 '^^ dfc' 

•^0 „.=-r^^^Ja Jo 



1 /-l p27ri(fe-fe')n 



n= — oo a, a' 



27ri 



= -TT-T / dec / dA; X 



C',fe(*) K',k{x)-^i^a,k e 



-27ri€Q. fct 



{x). 



The last line follows here from Poisson summation over n. 
In this expression d/dfc acts on a product of three terms, 
but as f ^ oo the dominant contribution comes from 
the derivative of the exponential. Neglecting the other 
two terms which are bounded, and using the orthonor- 
malization of Floquet states we finally obtain Eqs. (56), 
(57). 

For higher moments of the spatial distribution the ar- 
gument can be repeated and an analogous result is ob- 
tained 



a J" 

+o(t™-i) . 



dk 



(C5) 



